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ABSTRACT 

The AM Canum Venaticorum stars are rare interacting white dwarf binaries, whose 
formation and evolution are still poorly known. The Sloan Digital Sky Survey provides, 
for the first time, a sample of 6 AM CVn stars (out of a total population of 18) that 
is sufficiently homogeneous that we can start to study the population in some detail. 

We use the Sloan sample to 'calibrate' theoretical population synthesis models for 
the space density of AM CVn stars. We consider optimistic and pessimistic models 
for different theoretical formation channels, which yield predictions for the local space 
density that are more than two orders of magnitude apart. When calibrated with the 
observations, all models give a local space density po = 1 — 3 x 10~ 6 pc~ 3 , which is 
lower than expected. 

We discuss the implications for the formation of AM CVn stars, and conclude that 
at least one of the dominant formation channels (the double-degenerate channel) has 
to be suppressed relative to the optimistic models. In the framework of the current 
models this suggests that the mass transfer between white dwarfs usually cannot 
be stabilized. We furthermore discuss evolutionary effects that have so far not been 
considered in population synthesis models, but which could be of influence for the 
observed population. We finish by remarking that, with our lower space density, the 
expected number of Galactic AM CVn stars resolvable by gravitational- wave detectors 
like LISA should be lowered from current estimates, to about 1,000 for a mission 
duration of one year. 
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1 INTRODUCTION 



The Sloan Digital Sky Survey (SDSS; lYork et all l2000h 
is the largest photometric and spectroscopic sky survey 
ever carried out, consisting of imaging in five broad- 
band filters (u, g, r, i, z) and follow-up spectroscopy of 
a large fraction of targets outside the stellar (main- 
sequenc e) locus. The first phase, SDS S-I (i.e. Data Re- 
lease 5, lAdelman-McCarthv et al.l 120071 ) contains photom- 
etry of 8000 square degrees concentrated on the northern 
Galactic cap, down to a limiting magnitude of about g = 22, 
and spectroscopy of about 70% of that area (we shall re- 
fer to these as the 'photometric database' and the 'spectro- 
scopic datab ase', and we wil l use th e dereddened magnitudes 
according to ISchlegel et al.l (|l998T) throughout this paper). 
The spectroscopic database contains about 1 million objects, 
which have been obtained in about 1700 pointings ('tiles') 
using a 640-fibre multi-object spectrograph. Although the 
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spectroscopic follow-up to the imaging is mainly tailored 
towards finding quasars and galaxies, Galactic sources off 
the main sequence also make their way into the spectro- 
scopic sample, serendipitously or sometimes targeted as 
white dwarfs (WDs) and hot standard stars. 

Among the serendipitous targets are AM Canum Ve- 
naticorum stars: ultracompact binaries that consist of a 
white dwarf accreting helium-rich matter from a (semi- 
degenerate donor star. They are most easily picked up 
as helium-emission-line objects in the SDSS spectroscopic 
database. To date, six such syst e ms have been found 
in the SDSS (iRoelofs et all 120051 ; lAnderson et al. I 120051 : 
iGroot et al.ll200lih . a significant increase in the number of 
systems - only 10 AM CVn stars were known prior to the 
SDSS, and anot her two were recently discovered as s uper- 
nova candidates (|Woudt et alj|2005l : IPrieto et al.l l2006). The 
spectroscopic sample of emission-line systems from the SDSS 
is the first that is sufficiently large, complete, and homoge- 
neous that a study of the AM CVn population can be at- 
tempted. Since it is still a small sample, the approach in this 
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paper shall be to combine this observed sa mple with theo- 
retical pred ictions for the AM CVn stars l|Nelemans et al.l 
l200ll . 2004). in order to 'calibrate' the model populations. 

This will hopefully give us more insight into the as 
yet uncertain formation and evolution of these binaries. In 
addition, the AM CVn stars are of particular interest for 
gravitational- wave astrophysics: they are the only known ob- 
jects that can be detected with space-borne gravitational- 
wave detectors such as the propose d Laser Interferom- 
eter Space Antenna (LISA) mission ([Stroeer et al 1 120051 : 
IStroeer fc Vecchioll2006l ; IRoelofs et al.ll2007i ). Depending on 
their space density, as many as ~1 1,000 AM CVn stars can 
be res olved throughout the Galaxy by a LISA-type instru- 
ment l|Nelemans et alj|2004h . 



2 CALIBRATING POPULATION SYNTHESIS 
MODELS FOR THE AM CVN STARS 

2.1 Formation channels 

Three formation channels have been proposed to con- 
tribute to the population of AM CVn stars. In the first 
of these, the 'WD channel', a close pair of detached 
white dwarfs is formed via two common-envelope phases, 
which is subsequently driven towards Roche-lobe over- 
flow by gravitational- wave emission (e.g.lNather et al.lll98ll ; 
iTutukov fc Yungelsonl Il996l ; INelemans et al.l l200ll) . After 
the onset of Roche- lobe overflow at an orbital period of a few 
minutes, the system evolves towards longer orbital periods 
as the degenerate donor star expands upon mass loss. 

In the second scenario, the 'He-star channel', a white 
dwarf accretes from an initially non-degenerate helium star, 
while the binary evolves to shorter orbital periods. The he- 
lium burning in the core of the donor star quenches as a re- 
sult of the mass loss and the system reaches an orbital period 
minimum around ~10 minutes, at which point the donor star 
has become semi-degenerate and the s ystem starts to evolve 
towa r ds longer orbital periods again iTutukov fc Fedoroval 
1 19891 ; llben fc Tutukovlll99ll ; INelemans et al.ll200lh . 

The t hird, 'evolved-CV ch a nnel' was suggested more re- 
cently by IPodsiadlowski et alj (|2003T > and involves a Cata- 
clysmic Variable (CV) in which the mass transfer starts by 
the time the donor star starts to evolve off the hydrogen 
main sequence. Depending on the level of hydrogen depletion 
that is reached, such a system may shift its orbital period 
minimum from about one hour (the normal orbital period 
minimum for degenerate hydrogen-rich matter) down to the 
orbital period minimum for helium-rich matter at ~10 min- 
utes, with most such systems reaching intermediate orbital 
period minima due to remaining fractions of hydrogen in 
their cores. 

The objects emerging from each of these channels are 
expected to look like helium-rich, ultracompact AM CVn 
stars that evolve towards longer orbital periods and lower 
mass transfer rates, driven by angular momentum losses due 
to the emission of gravitational waves. 



2.2 Recent observational &: theoretical results 



Bilds ten et al.l l|2006l ) have recently shown that, at long or- 
bital periods, the AM CVn stars are expected to be dom- 



inated in the optical by an accretor that cools as if it 
were a stand-alone white dwarf. This is due to the rapidly- 
decreasing mass accretion rate losing control over the ther- 
mal state of the accretor, and the accretion luminosity 
dropping below the luminosity of the accretor. The steep 
decrease in mass transfer rate with orbital period, which 
is a consequence of the gravitational-wave losses being a 
steep function of orbital period, should cause the tran- 
sition from disc-dominated to accretor-dominated spectra 
to occur over a fairly small range in orbital period. Ob- 
servationally, this transition from a disc-dominated optical 
spectrum to an accretor-dominated spectrum is seen to oc- 
cur at orbital periods of about 30 minutes, with systems 
at longer orbital periods spending most of their time in 
a 'low' state characterized by helium emission lines from 
the disc and a continuum that is apparently dominated 
by the photosphere of the accretor. Systems at shorter or- 
bital periods spend most of their time in a 'high' state 
characterized by an absorption- line spectrum dominated by 
the accretion disc. T wo recent additions to the AM CVn 
family, SPSS J1240 (IRoelofs et all 120051') and 'SN2003aw' 
JWoudt fc Warned 120031 ; iNogami et al.ll2004l ; IRoelofs et ail 
l2006al ) at P — 37 and 34 minutes respectively, appear to 
spend most of their time in an emission-line state dominated 
by the accretor's photosphere, going into outburst not more 
than about once a year. Shorter-period systems such as CR 
Boo and V803 Cen at P — 25 and 27 minutes, on the other 
hand, spend most of their time in a high(er) state charac- 
terized by an accretion-disc-dominated absorption-line spec- 
trum; they are seen to spend less th an 10 per cent of their 
time in a low, emission-line state l|Patterson et al.1 Il997l . 
120001) . 

Assuming that, statistically, there is a transition from 
disc-dominated to accretor-dominated spectra at P m i n = 30 
minutes, and knowing the thermal evolution of the accretor- 
dominated systems, we can parametrize the colours and 
absolute magnitudes of these 'long-period' AM CVns at 
P > 30min with their orbital period. This allows us to de- 
termine the completeness of spectroscopic follow-up in the 
SDSS as a function of orbital period, by measuring the com- 
pleteness as a function of colour. If we further assume a cer- 
tain orbital period distribution for the AM CVn stars based 
on population synthesis studies, we can convolve this distri- 
bution with the completeness as a function of orbital period 
to obtain the total number of emission-line AM CVn stars 
present in the SDSS-I, from which we can derive the surface 
and space densities of these stars. 

Population synthesis studies of AM CVn stars have 
yielded predictions about their space density, but these 
are thought to be uncertai n to two orders of magnitude 
l|Nelemans et al.1 [200ll . 120041 ). Interestingly, the distribution 
of orbital periods is relatively well-known; it is mainly the 
absolute number of systems that is uncertain. One of the im- 
portant reasons for that is the fact that the angular momen- 
tum loss due to gravitational waves, which directly drives the 
evolution towards longer orbital periods in these systems, is 
a steep function of the orbital period. This causes the frac- 
tion of systems at long orbital periods (P > 30min), for 
instance, to be very similar for different population synthe- 
sis models since the evolution from the start of Roche-lobe 
overflow to this long-period regime is relatively instanta- 
neous - it takes only about 10 8 years. For all model popu- 
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lations as described in iNelemans et all (|200ll , l2004l ) we find 
only a ^ 2% fraction of systems at short orbital periods 
(P < 30min). 

This allows us to calibrate the total space density of 
AM CVn stars based on the emission-line population in the 
SDSS and the predicted orbital period distributions from 
population synthesis studies. To assess the influence of the 
exact shape of the orbital period distribution, we can calcu- 
late the total space density for different model populations 
and see how they differ. 

2.3 Formal method 

We wish to determine the space density of emission-line AM 
CVn stars from the observed N apec — 6 systems in the spec- 
troscopic database of SDSS-I. To this end we simulate ob- 
servations of model populations to see how many AM CVn 
stars we would have found in the SDSS-I if this population 
corresponded to the true population. The expected number 
of AM CVns in the spectroscopic database is given by an in- 
tegral over the sampled volume, multiplied by the modelled 
space density p' and by the completeness of spectroscopic 
follow-up C. The efficiency for identifying emission-line AM 
CVn stars once they have made their way into Sloan's spec- 
troscopic database is assumed to be (nearly) 100%: their 
characteristic helium-emission-line spectra are easily picked 
up using a combination of automated searches and visual 
inspection of spectra, as done by us and independently by 
lAnderson et~ai1 (|2005l ). The completeness of spectroscopic 
follow-up, i.e. the fraction of AM CVn stars in the photo- 
metric database of which a spectrum was taken, is expected 
to be far from 100%, and in addition it is expected to be a 
function of colour and apparent magnitude. One of our main 
concerns shall therefore be to determine this completeness. 

We express the expected number iV s 'p ec of emission-line 
AM CVn stars in the spectroscopic database of the SDSS, 
for a given model population, as an integral over the or- 
bital period P and the distance d, while summing over all 
spectroscopic pointings (primed quantities indicate model 
predictions) : 



oo d max (P) 




(1) 

Here, C(P, 6tiie, d) is the completeness of spectroscopic 
follow-up as a function of orbital period P (which we link to 
effective temperature and thus colour via a certain prescrip- 
tion), as a function of distance d (linked to apparent magni- 
tude g via a certain prescription) , and as a function of Galac- 
tic latitude b. We determine this completeness a posteriori 
from the Sloan databases by counting the fractions of ob- 
jects that have spectroscopic follow-up in colour, magnitude 
and Galactic latitude intervals. The quantity p' (P,b t n a ,d) 
is the population synthesis prediction for the space den- 
sity of systems as a function of orbital period, Galactic lat- 
itude and distance. The distance d max (P) is the distance at 
which an object of orbital period P would have an apparent 
magnitude <? max , where <? max is a suitably chosen magnitude 
limit. Finally, f2 t iio is the solid angle covered by each spec- 
troscopic tile (pointing), over which we take the sum to get 
the total area covered by the spectroscopic database. Thus 



the integral jQ,ti\ e d?Ad expresses the volume corresponding 
to a pointing. The solid angle covered by each pointing is 
slightly variable, but we will simply use an effective solid 
angle fitiic = f2 S pec/A?tiies, where O spcc = 5700 square de- 
grees is the total solid angle covered by the spectroscopic 
database, and iVtiies = 1700 is the number of spectroscopic 
tiles (pointings) over which we take the sum. 

The observed local space density po at Galactic height 
2 = 0, based on the observed N apec = 6, is then given by 



where p' is the local space density of the model population. 

The total number of emission-line AM CVn stars in 
the photometric database of SDSS-I down to g = (? max is 
obtained by setting C(P, fetiie, d) = 1 in eq.Q]and multiplying 
by f2 p hot/f2spec, where f2 p hot = 8000 square degrees is the 
solid angle on the celestial sphere covered by the photometric 
database: 



oo d max (P) 




This value N p hot divided by the area covered by the 
photometric database Q p hot gives the surface density on the 
sky a of emission-line AM CVn stars down to g — g max at 
the typical Galactic latitudes covered by the SDSS-I. 



2.4 Population synthesis models 

The double-degenerate WD channel and the single- 
degenerate He-star channel are thought to be the dominant 
contributors to the AM CVn population, with the evolved- 
CV c hannel contributing at the < 2% level l|Nelemans et al.l 
|2004 ). As in INelemans et all (|200lh we model two scenar- 
ios for the total population of AM CVn stars, an 'opti- 
mistic' and a 'pessimistic' model, each taking into account 
the two dominant channels. In the optimistic model, there 
is a strong tidal coupling between white dwarf binaries in 
the WD channel, which helps to stabilize the mass transfer 
by feeding back angular momentum accreted by the primary 
star to the orbit. Furthermore, carbon-oxygen (CO) white 
dwarfs accreting from a He star are not easily destroyed by 
edge- lit detonations (ELDs) in this model. In the pessimistic 
model we assume no tidal coupling between white dwarfs, 
causing many to merge and few to survive as stably accret- 
ing AM CVn stars. In this model we furthermore include a 
favourable scenario for ELDs to destroy AM CVn progeni- 
tors in the He-star channel, although this has only a limited 
effect on the number of AM CVn stars e merging from this 
channel. See INelemans et~ail |200ll . l2004h for a more quan- 
titative description of the two scenarios. 

In order to allow for the possibility that one of these 
two formation channels is shut off completely for some (as 
yet unknown) reason, we furthermore model optimistic and 
pessimistic populations for each of these channels individu- 
ally. 

All these population synthesis models include a mor e 
or less realistic Galactic model as in INelemans et al.l i|2004) . 
which assumes a star formatio n rate that depends on po - 
sition and time as modelled bv lBoissier fc Prantzosl |l999). 



4 G.H.A. Roelofs et al. 



All the stars that are formed are evolved using the latest ver- 
sion of the SeBa evolution code l|Portegies Zwart fc Verbunj 
1 19961 : iNelemans et al.ll200ll . 12004 ) . The resulting population 
of AM CVn stars at a distance R = 8.5 kpc from the Galactic 
centre and at time t = 13.5 Gyr is analysed. The only change 
we make to the model is that we assume a Galactic scale 
height of 300 pc for the AM CVn s tars as in Roelofs et al.l 
(2007), instead of the 200 pc used in lNelemans et all I^OoU T 
This is more appropriate for a population that is, on average, 
several Gyr old. 



2.5 Practical considerations 

The aforementioned parametrization of effective tempera- 
ture with orbital period P, from which we derive the colours 
as a function of orbital period, is taken to be 

T(P) = 19, 000 K x c -o.°25(p-p mm ) ^ 

based on figure 2 of iBildsten et all (|2006T ). with P in min- 
utes and P m in = 30min as before. From the same figure we 
parametrize the absolute magnitude as a function of orbital 
period as 



M g (P) = 10.5 + 0.075 (P - P m 



(5) 



with an estimated root-mean-square scatter of 0.5 magni- 
tude. There are two emission-line AM CVn stars for which 
there is a parallax measurement: GP Com and V396 Hya. 
GP Com at P = 46 m in has absolute magnitude My = 11.5 
i|Roelofs et al.l 120071 ) and V396 Hya at P = 65 min has 
My = 13.3 (J. Thorstensen, as quoted by iRoelofs et al.l 
2007). Both match well with the above parametrization. 

For the completeness C as well as for the model pop- 
ulation p' we have to construct a three-dimensional grid of 
values as a function of P, b and d. Since the modelled space 
density p' will (to a good approximation) only be a func- 
tion of P and Galactic height z = dsin(&), this reduces to 
a two-dimensional grid in P and z. We bin the modelled 
systems with a resolution of 5 minutes in P and 50 pc in z, 
from which we then fetch a population density at arbitrary P 
and z via bilinear interpolation. To give an idea of the model 
population, figure [T] shows the number of systems as a func- 
tion of orbital period P and Galactic latitude b, integrated 
up to (? m ax = 21. At high Galactic latitudes, the popula- 
tion is suppressed relative to the low-latitude population as 
an effect of the population scale height. This suppression is 
stronger at short orbital periods due to the systems' higher 
brightness at short periods, shifting the peak in the distri- 
bution to longer-period, fainter (and thus closer) systems. 
We see that above b ~ 30° , which is where the sky coverage 
of the SDSS-I is concentrated, the effect of Galactic latitude 
on the expected population is modest. 

The completeness of spectroscopic follow-up C(P,b,d) 
turns out to be insensitive to the Galactic latitude b over 
the area covered by the SDSS-I and we thus determine the 
completeness only as a function of colour, u — g and g — r, 
and apparent magnitude, g. Colour and apparent magnitude 
are again linked to orbital period and distance according to 
eqs.[4]and[5] See figure [2] for the completeness as a function 
of colour, integrated up to <? ma x = 21. To the lower left of 
the blackbody track in this figure we see the targeted DA 
white dwarfs, which deviate from blackbodies due to the 
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Figure 1. Modelled surface density tj' (per square degree) of 
AM CVn stars as a function of orbital period, down to g = 21, 
for three Galactic latitude ranges. The solid line shows the 'op- 
timistic' model, the dotted line the corresponding 'pessimistic' 
model. 
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Figure 2. Completeness of spectroscopic follow-up (grey-scale) 
as a function of u — g and g — r, to a limiting magnitude g = 21. 
Stars indicate the known AM CVns from SDSS-I. The solid line 
is the modelled cooling track for the AM CVns, closely following 
the blackbody cooling track (dotted line). Black dots indicate the 
'effective temperature' as goes into eq. [4] 



hydrogen absorption in their atmospheres, while the spec- 
troscopic follow-up to the right of the blackbody track is 
mostly targeted at quasars. The six AM CVn stars from the 
SDSS-I are seen to lie close to the blackbody curve; we model 
their cooling track parallel to but very slightly offset from 
it, so that it runs through the middle of the bunch of ob- 
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Figure 3. Completeness of spectroscopic follow-up C(P,g) as a 
function of the modelled orbital period P and for different inter- 
vals in apparent magnitude g. 

served systems. Figure [3] shows the completeness along the 
modelled cooling track as a function of apparent magnitude 
g. It is seen to depend strongly on apparent magnitude, and 
on colour in the 'cool' blackbody regime where the colours 
approach the bulk of the main sequence at P ~ 70 min. We 
construct a two-dimensional grid for the completeness C in 
intervals of 0.5 magnitude in g and 5 minutes in orbital pe- 
riod P from which we again get C at arbitrary P and g 
via bilinear interpolation. Since the SDSS has five colour 
bands, ugriz, the completeness of spectroscopic follow-up is 
not uniquely defined for each point in {u — g,g — r) space, 
but may in principle depend on the iz colours. In the hot 
blackbody regime where the AM CVn stars are, however, 
it is unlikely that the iz colours are of much influence. We 
verify this by determining the completeness as a function of 
r — i and i — z for the objects in a strip along the modelled 
cooling track from about 20,000-10,000 K, which includes 
the 6 AM CVn stars found in the SDSS. It turns out to be a 
very flat distribution; in other words the completeness along 
the cooling track in (u — g,g — r) is only weakly dependent 
on the (r — i, i — z) colours of the corresponding objects. 

Finally, we need to determine a suitable value g max to 
delimit a well-defined sample of AM CVn stars. We choose 
<?max = 21 for two reasons: first, the completeness of spec- 
troscopic follow-up in the SDSS plummets between g = 20.5 
and g = 21 (see also figure|3}, and second, up to this bright- 
ness limit we can still be confident that our efficiency for 
identifying emission-line AM CVn stars in the SDSS spec- 
troscopic database is (nearly) 100%. All six AM CVn stars 
identified in the spectroscopic database have g < 21. 

2.6 Results 

Table [T] shows the expected number of emission-line AM 
CVn stars in the SDSS spectroscopic database based on the 



different population synthesis models. From this we obtain 
the 'observed' local space density po by multiplying the mod- 
elled po with N Bpec /N Bpec (eq. 0, where N apcc = 6. We see 
that all models predict more systems than are observed, ex- 
cept for the pessimistic model in which additionally there 
are no AM CVns from the He-star channel. Although the 
space densities in the different models vary by more than 
two orders of magnitude, the observed space density based 
on these input populations varies only by a factor of 3, from 
po = 1 — 3 x 10 _6 pc -3 . As mentioned in section T2.2I this 
is because the variation in the total number of systems pro- 
duced in the different models is large whereas the variation 
in the orbital period distribution is quite small, and the lat- 
ter is what causes a variation in the calibrated (observed) 
space density. 

The uncertainties in the values in table[T]are dominated 
by the Poisson statistics of the observed sample (N apcc — 6) 
and the uncertainties in the parametrization of the tempera- 
ture and absolute magnitude with orbital period (eqs.|3]and 
[5]). The small sample of 6 systems contributes an intrinsic 
uncertainty of order 6 -1 ^ 2 tn 40% to the measured values in 
table [T] We shift the parametrization of absolute magnitude 
with orbital period (eq. [5j up and down by 0.3 magnitudes 
as an estimated uncertainty; this yields a variation of 32% in 
the observed space densities. Similarly we vary the effective 
temperature as a function of orbital period with ±20% and 
find a variation of 3-10% in the resulting space densities. 
Adding up all the numbers, we thus estimate a total uncer- 
tainty in the derived space densities that is close to a factor 

of 2 for each of the model populations. 

Note that the pessimistic model from iNelemans et al.l 

l|200lj ). when calibrated with the observed systems from the 
SDSS, actually yields a larger space density than the opti- 
mistic model. This is due to the fact that the pessimistic 
model is dominated by AM CVns from the He-star channel, 
which is on average a 'faster' channel that follows the star 
formation more closely. These systems are thus on average 
older and at longer orbital periods (and therefore fainter), 
causing the average completeness of spectroscopic follow-up 
to be lower and the resulting space density to be higher. 



3 DISCUSSION 

3.1 Previous space density estimates 

IWarnerl (|l995l ) estimates a local space density po ~ 2 x 
10 -6 pc -3 , identical to our result. Interestingly though, 
Warner arrives at this value for the wrong reasons: his esti- 
mate is essentially based on the short-period systems, and 
fails to take into account the large population of long-period 
systems that should accompany them (as mentioned before, 
the fraction of short-period systems is expected to be ^ 2%). 
But this is exactly offset by the assumption that the short- 
period systems have absolute magnitudes Mv = 9.5, which 
from HST parallax measurem ents we now know is actually 
M v ~ 6 (|Roelofs et alj|2007h . 

The estimate for the space density based on the short- 
period AM CVn stars was recently updated in lRoelofs et al.l 
<|2007f ) based on the updated absolute magnitudes from par- 
allax measurements, giving po ~ 1 x 10~ 6 /q pc -3 , where q 
is the known fraction of short-period systems with apparent 
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Model Modelled number Total in SDSS-I Modelled p' Observed po Observed Eq Per solar mass 





C-^spcc) 


(Afphot) 


(pc" 3 ) 


(PC" 3 ) 


(PC" 2 ) 


(^e 1 ) 


Optimistic 


106 


53 


2.7 x 10~ 5 


1.5 x 10~ 6 


4.5 x 10~ 4 


2.2 x 10- s 


Pessimistic 


12 


66 


6.1 x 10" 6 


3.1 x 10~ 6 


9.4 x 10" 4 


4.5 x 10- s 


He star only, opt. 


1(5 


68 


9.0 x 10" 6 


3.4 x 10~ 6 


1.0 x 10" 3 


4.9 x 10~ s 


He star only, pess. 


11 


67 


5.9 x 10" 6 


3.3 x 10" 6 


9.8 x 10" 4 


4.7 x 10~ s 


WD only, opt. 


90 


51 


1.8 x 10" 5 


1.2 x 10" 6 


3.5 x 10" 4 


1.7 x 10~ 5 


WD only, pess. 


0.82 


55 


2.2 x 10" 7 


1.6 x 10~ 6 


4.6 x 10" 4 


2.3 x 10~ s 



Table 1. O bserved space densities of AM CVn stars for differen t assumptions regardin g their populations. 'Optimistic' and 'pessimistic' 
models from INelemans etT al. (2001) with the Galactic model of iNelemans et all l l2004h . The observed local space density po is obtained 
by multiplying the modelled p' by N Bpcc /N!. pcc where N Bpcc = 6. The total A'phot is the number of emission-line AM CVn stars that 
should be present in the SDSS-I photometry down to g max = 21. The measured local surface density Sq projected onto the Galactic 
plane, in pc — 2 , provides a scale-height-independent space density. The number of AM CVns per unit mass (in Mg ) normalizes this 
number to the present-day, local stellar mass in the Galactic disk jChabrierll2"ooih . The observed numbers of AM CVns per pc and per 
Mg are accurate to an estimated factor of 2. The observed po in addition scales inversely with the assumed scale height; the numbers 
given are for a scale height of 300 pc. 



magnitude mv < 14.5. The space density derived here sug- 
gests that the unknown fraction of bright, short-period sys- 
tems is rather small. Alternatively, the comparable numbers 
could point towards a destruction of systems during their 
evolution from short to long orbital periods if one assumes 
that there should in fact be a significant population of short- 
period systems that have not yet been discovered. After all, 
the known AM CVns have been discovered serendipitously 
and predominantly at high Galactic latitudes, while the ma- 
jority of short-period systems with mv < 14.5 should in fact 
reside at low latitude, given their Mv — 6. 

3.2 Population synthesis results 

When comparing the observed space density to the models, 
one has to keep in mind that there are two different factors of 
uncertainty in the models: (a) the formation and evolution of 
AM CVn stars, and (b) the Galactic model containing the 
position- and time-dependent star formation history. This 
makes it more difficult to draw conclusions about one of the 
two based on the observed space densities. 

Fortunately the uncertainties due to the Galactic model 
are expected to be limited. We are looking at the local pop- 
ulation of AM CVn stars, and what matters thus is the lo- 
cal star formation history, which is calibrated observation- 
ally with for instance the local stella r density, giving the 
integ rated star formation history (see IBoissier fc Prantzosl 
1 19991 ) . Inclusion of a star formation history that is a function 
of both time and distance from the Galactic centre, versus 
one where the time-dependence is independent of position, 
has lowere d the local space dens i ty of AM CVn stars by a 
factor 2.7 (|Nelemans et al.ll200ll . |2004| ). This is due to the 
fact that the AM CVn star progeny are relatively old popula- 
tions; the delay in the star formation at R — 8.5 kpc relative 
to the centre of the Galaxy reduces the number of systems 
in the solar neighbourhood that have already evolved to be- 
come AM CVn stars. We estimate that the remaining room 
for variations in the modelled space density of AM CVn stars 
due to the local star formation history is less than a factor 
of 2. 

What could further be of influence is the assumed scale 
height h z = 300 pc of the AM CVn stars above the Galactic 



plane. Figure [4] shows the modelled and observed distribu- 
tions in apparent magnitude g for different scale heights of 
h z = 200, 300 and 400 pc. The effect on the total number of 
expected systems in the SDSS-I is clearly very small, indi- 
cating that the discrepancy between the observed and the 
modelled space densities is not due to the assumed scale 
height of the model population. Note also that the larger 
scale heights appear to be slightly (though not significantly) 
preferred based on the observed sample. Although the ratio 
between the modelled and the observed space densities is 
insensitive to the assumed scale height, the absolute num- 
bers are: if assuming a scale height of 400 pc instead of 300 
pc, both numbers need to be multiplied by approximately 
3/4. In order to provide a scale-height-independent number, 
which is essentially what we obtain from the SDSS-I since 
we are looking out to several scale heights for the average 
AM CVn star, we provide in table Q] the surface density £ 
of AM CVns when we project them all onto the Galactic 
plane, in units of pc -2 . Adopting a local surface density of 
21Ma /pc" 2 for stellar matter in the Galactic disk l|Chabrier| 
2001), we also give the local number of AM CVns per solar 
mass, in the last column of table [T] 

We conclude that a substantial part of the overopti- 
mism in the 'optimistic' model (which yields a factor 18 
higher space density than observed, see table[TJ is due to the 
evolutionary model of the AM CVn stars. The easiest expla- 
nation is that the optimistic model for the WD channel is 
far too optimistic: mass transfer between white dwarfs upon 
Roche-lobe overflow cannot be stabilized in the majority of 
cases , but rather leads to merger events (INeleman s et al.l 
l200ll : lMarsh et al.ll20oj ). The pessimistic model for the WD 
channel however, which assumes no stabilizing effect on the 
mass transfer from feedback of angular momentum due to 
tidal coupling between white dwarfs, yields predictions for 
the space density that are too low. This means that there 
is either a contribution from the He-star channel, which by 
itself produces modelled space densities that are compatible 
with the observed value, or that there may be some stabi- 
lizing effect due to tidal coupling in the WD channel. 

The two most viable scenarios in the framework of the 
current models are therefore: 
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Figure 4. Modelled distribution of spectroscopically identified 
AM CVn stars in SDSS-I as a function of apparent magnitude 
g. The dotted, solid and dashed lines represent model population 
scale heights of 200, 300 and 400 pc, respectively. The upper set 
of lines are for the optimistic model; the pessimistic model (lower 
lines) gives virtually identical results but scaled downwards. 



(i) There is a mild tidal coupling between white dwarfs, 
which helps to stabilize the mass transfer in up to 10% of 
all cases in which this could work, while the He-star channel 
is inactive; 

(ii) The He-star channel is active and there is no efficient 
tidal coupling between white dwarfs in the WD channel, 
causing the WD channel to be of little influence. 



3.3 Implications for AM CVn star evolution 

If we look outside the box of the current models, the rela- 
tively low number of observed emission-line AM CVn stars 
may most readily be explained by assuming that long-period 
AM CVn stars do not always exhibit the helium emis- 
sion lines by which they were identified from the SDSS-I. 
It has been sugges ted from donor star model calculations 
l|Delove et al.l 120051 ) that in AM CVn stars with hot, semi- 
degenerate donors, these donors may cool and contract, so 
that the mass transfer rate goes down substantially and, de- 
pending on how the timescales work out, they might even be- 
come detached. This would eliminate them from the sample 
studied here. Essentially, such semi-degener ate donors cool 
down when their (large) therma l timescale (|Savoniie et al.l 
ll98rj ; |Tutukov fc Fedorowalll989l ) finally catches up with the 
(rapidly increasing) timescale of angular momentum loss 
from the system due to the emission of gravitational waves. 
Calculations of the thermal evolution of the donor stars in 
AM CVns (Deloye in preparation; private communication) 
indicate that this could happen at P « 45 min. Observa- 
tionally, it has recently been found that the donor stars in 
short-period AM CVns (P < 30 min) appear to be mostly 
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Figure 5. Modelled distribution of spectroscopically identified 
AM CVn stars in SDSS-I as a function of orbital period, along 
our modelled cooling track. The 4 Sloan AM CVns for which a 
measured orbital period exists are indicated in the figure. The 
solid and dotted lines are the optimistic and pessimistic models, 
respectively. 



semi-degenerate l)Roelofs et al ] |2007l : lNasser et alJlioOlh ; it 
would thus appear that this scenario is possible. 

Observationally, this would become clear as a fall-off 
of the population at P > 45 min or a pile-up of systems 
around P ~ 45 min, depending on whether they become 
fully detached. The current sample of six, of which orbital 
periods have only been measured for four, at present pre- 
vents us from drawing conclusions. The fact that AM CVns 
at orbital periods well above 45 minutes have been observed, 
indicates that at least not all of them become detached in- 
definitely. In addition, even if this were the case, it is clear 
from figure[2]that this could at most have caused a suppres- 
sion of the population by a factor of 3, based on the fraction 
of observed systems that is expected to be at orbital periods 
P > 45 min. 

The He-star models take into account optimistic and 
pessimistic estimates for the destruction of AM CVn stars 
by edge-lit detonation s of the accreting CO white dwarf 
l|Nelemans et al.l 1200 ll ). Although recent studies question 
the e xistence of such ELDs altogether (e.g. lYoon fc Langerl 
|2004| ). their possible non-existence is not very relevant for 
the AM CVn population since it is of relatively little in- 
fluence on the number of AM CVn stars produced in the 
He-star channel (see table [TJ. What could be of influence is 
the suggestion that the He-star channel may hardly work at 
all due to the poor efficiency for ejecting the common enve- 
lope in such a configuration. This is thought to be due to the 
hydrodynamics of the relatively weak density gradient be- 
tween the He-burning core and its hydroge n mantle, which 
togeth er form the He-star progenitor (see Sandquis t et al.l 
(2000) for a study of how this may prevent the formation of 
low-mass (< 0.2 Mq) He WDs from a common envelope). 



8 G.H.A. Roelofs et al. 



In a variation on the ELD theme, iBildsten et al.l (|2007f ) 
have recently shown that a substantial fraction of AM CVns, 
and possibly all of them, could host a strong thermonuclear 
event due to the ignition of a sufficiently thick layer of ac- 
creted helium on the accretor's CO core. This could destroy 
the binary and thus remove AM CVns from the population 
after they have started their evolution towards longer orbital 
periods as stable mass-transferring binaries. 

In the WD channel, the effect of the unknown tidal cou- 
pling between the accretor and the donor on their survival 
as stable mass-transferr ing white dwarfs jNelemans et alj 
l200ll ; iMarsh et al.l 120041 ') has been taken into account by 
modelling an optimistic (perfect coupling) and a pessimistic 
(no coupling) scenario. In the latter case, the WD chan- 
nel is so much suppressed due to the mass transfer becom- 
ing unstable that the observed space density is too high for 
this channel to be the only contributor. An additional, rel- 
atively unstudied effect is the possible ignition of He-core 
white dwarfs upon accretion of a substantial amount of he- 
lium from their helium- white-dwarf do nors, either through a 
series of inward-moving helium flashes i|Saio fc JeffenfeOOOT i 
or via an inward-moving burning front (Bildsten in prepa- 
ration; private communication). Such He+He WD binaries 
represent up to 20% of the AM CVn star progenitors in the 
WD channel for the optimisti c mod el (but note again the 
arguments bv lSandauist et alj (2000) against the formation 
of low-mass He WDs in such binaries); the possible non- 
existence of AM CVns descending from He+He WD binaries 
should therefore be of limited influence. In the pessimistic 
model, only WD binaries of extreme mass ratio survive the 
initial mass transfer, leaving no AM CVn stars from He+He 
white dwarf binaries anyway. 

With the lower space densities reported in this pa- 
per, we note that the e volved-CV channel for the forma- 
tion of AM CVn stars (|Podsiadlowski et all l2003h . which 
was previously estima ted to contribute at the < 2% level 
l|Nelemans et al.ll2004T ). might become important. However, 
AM CVns from this channel should either (a) pile up more 
strongly at long orbital periods due to the fact that most 
systems are not sufficiently hydrogen-depleted to shift their 
orbital period minima all the way down to ~10 minutes; or 
(b) look like 'ordinary' CVs due to their remaining hydro- 
gen, implying that there should be many CVs with orbital 
periods well below one hour. Neither of that appears to be 
observed. Nevertheless the observational challenge of finding 
an AM CVn star with traces of hydrogen in its accretion disc 
remains. 

3.4 Implications for the AM CVn stars as LISA 
sources 

AM CVn stars are thought to be among the most important 
sources for future space-borne gravitatio nal-wave detectors 
such as LISA (e.g. iNelemans et al.ll200"3 ). In addition, it is 
currently the only c lass of sources for which known, resolv- 
able sources exist (|Stroeer et all 120051 ; IStroeer fc Vecchiol 
l200d ; lRoelofs et al.ll2007l ). Current estimates for the number 
of AM CVn stars resolvable with a LISA-type gravitational- 
wave detector are based o n the 'optimistic ' mode l for their 
formation and evolution (|Nelemans et al.l [2004). The re- 
sults presented here suggest that these estimates should 
be lowered. As discussed in the previous section, a loss 



of systems during their evolution from the most sensitive 
LISA band (P < 20min) to the emission-line regime at 
P > 30 min could cause us to underestimate the short- 
period population. Nevertheless it would seem reasonable, 
based on the results shown in table [TJ to lower the estimate 
of ~ 1 1,000 Galactic AM CVn stars that can be resolved 
with LISA to about ^1,00 0, for a mission duration of one 
year l|Nelemans et al.ll2004T ). 



4 CONCLUSION 

We have investigated the population of helium-emission- 
line AM CVn stars based on the sample of six new sys- 
tems from the Sloan Digital Sky Survey. We have com- 
pared the observed sample to predicted samples as obtained 
from population synthesis models. For different model pop- 
ulations we have derived observed local space densities of 
1 — 3 x 10 -6 pc -3 . In addition to the variations in observed 
space density for the different models, we estimate uncer- 
tainties of a factor of 2 due to Poisson noise and (assumed) 
uncertainties in the colours and absolute magnitudes of the 
emission-line AM CVn stars. In the standard evolutionary 
picture, where all AM CVn stars evolve to become emission- 
line systems at orbital periods above ~30 minutes, this cor- 
responds to a space density of the entire population that 
is very nearly identical to this value (within 2%). A fur- 
ther refinement of the space density would require addi- 
tional follow-up spectroscopy of photometric candidates in 
the SDSS, to increase the spectroscopic sample and bring 
down the Poisson noise. Since there are relatively few other 
types of objects with approximate hot blackbody colours 
(mainly single DB white dwarfs), this is feasible. Further 
orbital period measurements would be needed to reduce the 
uncertainty in the colour evolution and orbital period dis- 
tribution of the AM CVn population. 

The observed space density presented here is lower than 
expected from the population synthesis models, which pre- 
dict 6 x 10 -6 pc -3 for the pessimistic models to 3 x 10 -5 pc -3 
for the optimistic models. At least one of the proposed dom- 
inant formation channels, the double-degenerate WD chan- 
nel, has to be suppressed by at least an order of magni- 
tude compared to the optimistic models, which assume that 
the mass transfer between white dwarfs can be stabilized 
in many cases due to a strong tidal coupling of spin and 
orbital angular momentum. A significant effect of tidal cou- 
pling on the survival rate of AM CVns from the WD channel 
is possible only if the second dominant formation channel, 
the single-degenerate He-star channel, is also severely sup- 
pressed. 

We have presented an inventory of ideas, based on cur- 
rent theory, of how the formation channels may be sup- 
pressed relative to the models we have used. Most of these 
(detachment of long-period AM CVns, ignition of He+He 
WD accretors) have only a limited effect, but inefficient 
common-envelope ejection in the He-star channel could po- 
tentially be effective in shutting down the He-star channel 
completely. 



The population of AM CVn stars 9 



5 ACKNOWLEDGMENTS 

It is a pleasure to thank Jim Liebert and Chris Deloye for 
stimulating discussions. GHAR and PJG were supported by 
NWO VIDI grant 639.042.201 to P.J. Groot. GN was sup- 
ported by NWO VENI grant 639.041.405 to G. Nelemans. 

Funding for the SDSS has been provided by the Al- 
fred P. Sloan Foundation, the Participating Institutions, the 
National Science Foundation, the U.S. Department of En- 
ergy, the National Aeronautics and Space Administration, 
the Japanese Monbukagakusho, the Max Planck Society, 
and the Higher Education Funding Council for England. The 
SDSS is managed by the Astrophysical Research Consortium 
for the Participating Institutions. The Participating Institu- 
tions are the American Museum of Natural History, Astro- 
physical Institute Potsdam, University of Basel, Cambridge 
University, Case Western Reserve University, University of 
Chicago, Drexel University, Fermilab, the Institute for Ad- 
vanced Study, the Japan Participation Group, Johns Hop- 
kins University, the Joint Institute for Nuclear Astrophysics, 
the Kavli Institute for Particle Astrophysics and Cosmol- 
ogy, the Korean Scientist Group, the Chinese Academy 
of Sciences (LAMOST), Los Alamos National Laboratory, 
the Max-Planck-Institute for Astronomy (MPIA) , the Max- 
Planck-Institute for Astrophysics (MPA) , New Mexico State 
University, Ohio State University, University of Pittsburgh, 
University of Portsmouth, Princeton University, the United 
States Naval Observatory, and the University of Washing- 
ton. 



References 

Adelman-McCarthy J.K., et al., 2007, ApJS, in press 

(|arXiv:0707.3380V l) 
Anderson S.F., et al., 2005, AJ, 130, 2230 
Bildsten L., Shen K.J., Weinberg N.N., Nelemans G., 2007, 

ApJL, 662, 95 

Bildsten L., Townsley D., Deloye C.J., Nelemans G., 2006, 

ApJ, 640, 466 
Boissier S., Prantzos N., 1999, MNRAS, 307, 857 
Chabrier O, 2001, ApJ, 554, 1274 

Deloye C.J., Bildsten L, Nelemans O, 2005, ApJ, 624, 934 
Groot P.J., Roelofs G.H.A., Woudt P., Nelemans G., 
Liebert J., Steeghs D., Marsh T.R., Warner B., 2006, MN- 
RAS, submitted 
Iben I. Jr., Tutukov A.V., 1991, ApJ, 370, 615 
Marsh T.R., Nelemans G., Steeghs D., 2004, MNRAS, 350, 
113 

Nasser M.R., Solheim J.-E., Semionoff D.A., 2001, A&A, 
373, 222 

Nather R.E., Robinson E.L., Stover R.J., 1981, ApJ, 224, 
269 

Nelemans G., Portegies Zwart S.F., Verbunt F., Yungelson 

L.R., 2001, A&A, 368, 939 
Nelemans G., Yungelson L.R., Portegies Zwart S.F., 2004, 

MNRAS, 349, 181 
Nogami D., Monard B., Retter A., Liu A., Uemura M., 

Ishioka R., Imada A., Kato T., 2004, PASJ, 56, 39 
Patterson J., et al., 1997, PASP, 109, 1100 
Patterson J., Walker S., Kemp J., O'Donoghue D., Bos M., 

Stubbings R., 2000, PASP, 112, 625 



Podsiadlowski Ph., Han Z., Rappaport S., 2003, MNRAS, 
340, 1214 

Portegies Zwart S.F., Verbunt F., 1996, A&A, 309, 179 
Prieto J., Anderson S., Becker A., Marriner J., Sako M., 

Jha S., 2006, Central Bureau Electronic Telegrams, 692, 

1 

Roelofs G.H.A., Groot P.J., Marsh T.R., Steeghs D., Barros 
S.C.C, Nelemans G., 2005, MNRAS, 361, 487 

Roelofs G.H.A., Groot P.J., Marsh T.R., Steeghs D., Nele- 
mans G., 2006a, MNRAS, 365, 1109 

Roelofs G.H.A., Groot P.J., Marsh T.R., Steeghs D., Nele- 
mans G., 2006b, MNRAS, 371, 1231 

Roelofs G.H.A., Groot P.J., Benedict G.F., McArthur B.E., 
Marsh T.R., Steeghs D., Morales-Rueda L., Nelemans G., 
2007, ApJ, 666, 1174 

Saio H., Jeffery C.S., 2000, MNRAS, 313, 671 

Sandquist E.L., Taam R.E., Burkert A., 2000, ApJ, 533, 
984 

Savonije G.J., de Kool M., van den Heuvel E.P.J. , 1986, 
A&A, 155, 51 

Schlegel D.J., Finkbeiner D.P., Davis M., 1998, ApJ, 500, 
525 

Stroeer A., Vecchio A., 2006, Class. Qu. Grav., 23, 809 
Stroeer A., Vecchio A., Nelemans G., 2005, ApJ, 633, 33 
Tutukov A.V., Fedorova A.V., 1989, Soviet Astronomy, 33, 
606 

Tutukov A.V., Yungelson L., 1996, MNRAS, 280, 1035 
Warner B., 1995, Cataclysmic Variable Stars. Cambridge 

Univ. Press, Cambridge 
Woudt PA., Warner B., 2003, MNRAS, 345, 1266 
Woudt PA., Warner B., Rykoff E., 2005, IAU Circular, 

8531, 3 

Yoon S.-C, Langer N., 2004, A&A, 419, 654 
York D.G., et al., 2000, AJ, 120, 1579 



